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DYNAMICAL UPPER BOUNDS ON WAVEPACKET SPREADING 


ROWAN KILLIP, ALEXANDER KISELEV, AND YORAM LAST 


Abstract. We derive a general upper bound on the spreading rate of wavepackets 
in the framework of Schrodinger time evolution. Our result consists of showing 
that a portion of the wavepacket cannot escape outside a ball whose size grows 
dynamically in time, where the rate of this growth is determined by properties 
of the spectral measure and by spatial properties of solutions of an associated 
time independent Schrodinger equation. We also derive a new lower bound on 
the spreading rate, which is strongly connected with our upper bound. We apply 
these new bounds to the Fibonacci Hamiltonian—the most studied one-dimensional 
model of quasicrystals. As a result, we obtain for this model upper and lower 
dynamical bounds establishing wavepacket spreading rates which are intermediate 
between ballistic transport and localization. The bounds have the same qualitative 
behavior in the limit of large coupling. 


1. Introduction 

Let R be a self adjoint operator on a separable Hilbert space 7Y. The time- 
dependent Schrodinger equation of quantum mechanics, = n't!), leads to a uni¬ 
tary dynamical evolution in Ti, given by 

( 1 . 1 ) m = 


We are primarily interested here in cases where R is a Schrodinger operator of the 
form A -|- U on a space such as or and where the initial vector is 

a localized wavepacket. Under the evolution ( |1 . 1| ) , the wavepacket will often spread 
with time and it is the nature of this spreading which is our main object of interest 
here. This question, “What determines the spreading of a wavepacket” , has been 
an active held of research over the last 15 years and there is by now a considerable 

0 |, |, 0 , 0 , 0 0 , 0 , 0 , 0 , 0 


body of literature devoted to it 





m m 


is just a partial list). We note that there are 
some situations, such as some systems with absolutely continuous spectrum which 
can be studied by scattering theory ||^ , where good understanding of the dynamics 


exists and is much older. The primary focus of the more recent studies has been 
on greater generality and, in particular, on systems exhibiting spectral phenomena 
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such as singular continuous spectrum and dense point spectrum. In particular, sin¬ 
gular continuous spectrum tends to occur for basic models of quasicrystals (see, e.g., 
||T^ , |i0[| ) and the study of wavepacket spreading in such models is important to the 
understanding of transport properties in such media. 

A major focus of many recent studies has been to elucidate the relations between 
wavepacket spreading rates and continuity properties of spectral measures. Recall 
that each initial vector -0 has a spectral measure /i^, dehned as the unique Borel 
measure obeying 'll:) = f{x) d^^{x) for any measurable function /. Here 

we use (■, ■) to denote the scalar product in H. A major discovery of Guarneri |^, p!7| , 
which has also been extended by several other authors B1.1.0 0.0. is that 
appropriately measured continuity properties of the spectral measure imply fairly 
general lower bounds on the spreading rate of the wavepacket. Recently, Kiselev-Last 
p5| (also see for a related heuristic result and for a recent extension) proved 
a variant of Guarneri’s bound which blends continuity properties of the spectral 
measure with the spatial decay rate of generalized eigenfunctions. This bound is 
generally stronger than what one can obtain from properties of the spectral measure 
alone and it is also of somewhat broader applicability. We note that while there 
is some pure theoretical interest in relating wavepacket dynamics to such things as 
continuity properties of spectral measures and spatial properties of eigenfunctions, the 
interest in the above bounds is more than purely theoretical. Properties of spectral 
measures and eigenfunctions can be studied for concrete models of interest and the 
above discussed bounds have been used to obtain dynamical results in cases where 
there is currently no alternative approach to study the dynamics. In particular, this 
approach has been used to obtain dynamical results for a number of quasiperiodic 
operators B00I and for operators with decaying potentials p5| . 

The situation with putting upper bounds on wavepacket spreading rates is much 
more problematic. There is a fairly general argument of Simon 


showing that for a 

broad class of Schrodinger operators (in particular, every operator of the form A + V 
on £^(Z'^)) wavepacket spreading rates cannot be more than ballistic (namely, linear 
in time). For the regime of dynamical localization (namely, situations where the 
wavepacket does not spread beyond a hnite region of space and so, in particular, the 
spectrum must be pure point) there is a fairly good understanding |T^, [T^, of how 
suitable spatial properties of eigenvectors (or of spectral projections) are related to 
dynamics and how to specify necessary and sufficient conditions for the occurrence of 
complete dynamical localization. For problems with continuous spectrum, however, 
it had been a noted open problem to get any results beyond Simon’s ballistic bound, 
namely, to be able to specify conditions that would guarantee wavepacket spreading 
rates that are slower than ballistic. 

The main aim of the current paper is to provide a certain solution to this problem. 
What we are able to do, in essence, is to specify conditions that would guarantee that 
some portion of the wavepacket does not spread too fast (namely, it remains within 
a box whose size grows with time at an appropriate rate). We cannot control the 
entire wavepacket and thus we cannot bound the growth rate of quantities that are 
sensitive to small, fast spreading portions of the wavepacket, such as moments of the 
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position operator. Nevertheless, we believe that our result is in many ways a natural 
complement to some of the Guarneri type lower bounds discussed above. We note that 
while Guarneri type bounds, such as the variants obtained in |]^, |?7||, have been often 


formulated in terms of moments of the position operator, such bounds are usually 
obtained by controlling only a portion of the wavepacket. It so happens that in order 
to establish a lower bound on the growth rate of moments of the position operator it 
suffices to show that some portion of the wavepacket is spreading at an appropriate 
rate, whereas bounding such quantities from above would require tight control of the 
entire wavepacket—a much more delicate task. Moreover, we will see that our upper 
bound has a strong connection with the Kiselev-Last |2^ lower bound discussed 


above. While we formulate things somewhat differently here, the basic information 
going into our upper bound is, roughly speaking, the same combination of spectral 
measure continuity and spatial behavior of generalized eigenfunctions which goes into 
that lower bound. Indeed, we also formulate here a new stronger variant of this lower 
bound showing that, in essence, we have good understanding of how the spreading 
rate of the slowest moving portion of the wavepacket is determined. The behavior of 
the fastest moving portion remains outside our scope of control and, in particular, 
establishing sub-ballistic upper bounds on the growth rate of moments of the position 
operator remains an interesting open problem. 

From here on we study tridiagonal operators of the form 


( 1 . 2 ) 


{Hu){n) = a{n)u{n + 1) + a{n — l)u{n — 1) + b{n)u{n) 


(also called Jacobi matrices) dehned on £^(Z_|_) or Here h[n) and a{n) are real 

numbers, and a{n) ^ 0 for any n. Moreover, we assume that 

which is sufficient to ensure that these operators are self-adjoint (see, e.g., 0). We 
note that discrete one-dimensional Schrodinger operators of the form Aj-K on £^(Z_|_) 
or £^(Z) are a special case of such tridiagonal operators (obtained by setting a{n) = 1 
and h{n) = V{n)). While such one-dimensional operators will be our main interest 
here, we discuss in Section 7 below how our results are also applicable to more general 
problems and, in particular, how our dynamical upper bound is also applicable to 
multidimensional Schrodinger operators of the form A -|- K on £^(2*^). 

We will look at the evolution of the initial vector = Ji (where Jn(fc) = Jnfc), but 
the results can also be recast for other (? vectors. 

Let ueiji, z) be the solution of the equation 

(1.3) a{n)u{n + 1) + a{n — l)u{n — 1) + h{n)u{n) = zu{n), 


z G C, satisfying ug(0, z) = sin 6*, ^^(l, z) = cos 9. Given a function 0 : Z_|_ —>• C, we 
dehne for any L > 0, 


ml 


LG 


n=l 


(1.4) 
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where we use the convention that the sum is zero if the upper limit is less than the 
lower one. Similarly, for a function 0 : Z — >■ C, we dehne 

lL2i 

= E + Vi - - i)t + (T - [L,m[L2i + i)L 

n=-[LiJ 

Here [a:J is the greatest integer less than or equal to x. While the extension to real L 
is not essential, it will be technically convenient. Consider first the half-line problem, 
namely, the operator on £^(Z+). Without loss of generality, we restrict consideration 
to a Dirichlet boundary condition {6 = 0), where the operator o is defined on 
functions vanishing for n < 0. We denote the resulting self-adjoint operator by 
We also assume here a(0) = 1. Note that the operator itself is independent of 
a(0). The choice of a(0) is only effecting the definition of the solutions uo{n, z) for 
6^0. Given e > 0 and an energy G M, we dehne the length scale L'^{E) by 

(1.5) max||M 0 (n,E)||^+(g) ■ min \\ue{n, = e“^ 


The left-hand side is monotonely increasing as the region of summation expands 
and so L'^iE) is well-dehned for every e > 0. We have L'^{E) —>• 1 as e —>• cx) and 
L^{E) —>■ cx) as e —> 0. Let m+{z) be the Weyl-Titchmarsh m function corresponding 
to the operator 

m+{z) = ((iL+ - z)~^6i, (5i) 


The spectral measure /r'*' of the operator is dehned by the equality 


m+{z) 


X — z 


and coincides with fiSi, the spectral measure of the vector hi. As in previous works 
on the subject, we deal with time-averaged quantities. For any function A{t) of time, 
we denote 


{A(t))T = ^ Je A{t) dt. 
0 


This type of averaging is slightly different from the Cesaro averaging used in many 
previous works and it is more convenient for what we do here. Note that the difference 
is not very signihcant. In particular, power law behaviors in T must be the same for 
both kinds of average. 

Our first result is the following 


Theorem 1.1. Let he the half-line operator defined on by (H) and a 

Dirichlet boundary condition. Let the characteristic scale L^_i (E) be defined by (0). 
Then for any T > 0 and L > 1, we have 

(1.6) (l|e-"+‘ii||i)T > ({£| L+_,(B) < L}) , 

where C is a universal positive constant. 
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Remarks. 1. The expression on the left hand side of ( |1.6| ) gives the averaged norm 
of the portion of the wavepacket remaining in a ball of size L. Theorem |1.1| shows 
that this norm is bounded from below by a quantity proportional to the norm of the 
spectral projection of (5i on the set of energies where L^_i{E) < L. Thus, if we choose 
the size L to be greater than a certain scale which depends on time and properties of 
the solutions uq, we are guaranteed to have a significant portion of the wavepacket 
remaining (on average) in a ball of this size. 

2. While the universal constant C can be explicitly estimated from our proof below, 
we made no real effort to obtain an optimal value for it. Our technique can only 
yield a number which is significantly smaller than 1 and thus (|1.6|) can control the 
location of only a portion of the wavepacket. 


The next criterion is a simple corollary of Theorem 1.1 


It relates an upper bound 
on wavepacket spreading to the growth of the norms of transfer matrices. Let us 
denote by <l>(n, E) the transfer matrix from the site 0 to n : 


<h(n, E) 


Mo(n + 1, E) M^/2(n + 1, A 
Uoin^E) M^/2(n,E) )' 


We also denote, for L > 1, 


LLJ-l 

n=l 


where ||$(ri, i?)|| is the usual operator norm of the matrix <l)(n, E). Define L^{E) by 
(1") ||4(i5)||y(,5, =2||4(1,B)-'||€-'. 

Then we have the same result for a new characteristic scale: 


Corollary 1.2. Let he the half-line operator defined on by (|1.2| ) and a 

Dirichlet boundary condition. Let the characteristic scale Llf,_fiE) he defined by (|1.7]). 
Then for any T > 0 and L >2, we have 

(1,8) (l|e-"+‘ii||!)T > ({E\Lp,(E) < i}) . 

where C is a universal positive constant. 


Remark. By the Landauer formula |^, the resistance p{n) of a sample of size n is 
given by 


p(n) = i[iTr($*(n,E)«h(n,E))-l] 

(where is the transpose of $). Hence, the growth rate (with L) of ||<h(i?)||| is 
connected with the growth rate of the sum 


L 

5 ^ pH- 

n=l 






6 


ROWAN KILLIP, ALEXANDER KISELEV, AND YORAM LAST 


Corollary |L^ can thus be viewed as a confirmation of a physically reasonable state¬ 
ment that higher resistance leads to slower transport. 

We will show below that L'^{E) < L+(£'), and hence Corollary |1.2| is strictly 


weaker than Theorem |1 . 1| . However, it may be more straightforward to apply it in 
some situations. 

In addition, we derive a new lower bound on wavepacket spreading which is also 
related to the characteristic length scale L'^{E). 

Theorem 1.3. Let he the half-line operator defined on by o with a 

Dirichlet boundary condition. Denote by Ps the spectral projection on some Borel set 
S. Then for every T > 0 and L > 0, 


(1.9) 


\e-^^^^PsSi\\l)T<C 


\\Mn,E)\\] 






where C is a universal constant. 


Roughly, Theorem 1.3 shows that if L is such that the ratio 


\\uo{n,E)\\l/\\uo{n,E)\\l, 

is small for E E S, then most of the portion of the wavepacket which is supported on 
energies within the set S leaves the ball of radius L by the time T. 

While we dehned the scales Lf (E) purely in terms of solutions of the equation ( p..3|) , 
we will see below that the quantities which enter in this definition also determine the 
local scaling properties of the spectral measure. This connection can be turned around 
to control Lf{E) using partial information on solution behavior along with continuity 
properties of the spectral measure. In particular, the meaning of Theorem and its 
connection with the Kiselev-Last lower bound can be clarified by noting that it 
implies the following. Let a{E) be a measurable function such that for a.e. E w.r.t. 


( 1 . 10 ) 


= limsup 

^ ^ (2e)“(®) 


< CX). 


Let '^{E) be a measurable function such that for a.e. E w.r.t. pf 


( 1 . 11 ) 


limsupL '^^^"^\\uQ{n^E)\\] 


< CX). 


We will show that Theorem |1.3| implies: 

Theorem 1.4. Leta{E), '^{E) satisfy (|1.10|) , (|1.11|) and let ri{E) = a{E)/'j{E). Let 
S C M. be a Borel set and assume that for all energies E E S, r}{E) > 6 > 0. Then 
for any (/ > 0, there exists a constant Cg such that for all T > 0, 

(1.12) (||e-^+^P55i|lLT^)T < 1 - /X+ {{E I p{E) > b}) + g. 


Remarks. 1. As usual |^, from the estimate ( |1.12|) follow lower bounds on the 
growth rate of moments of the position operator. 

2. Theorem |1.4| is a local version (and thus also a somewhat stronger variant) of 
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Theorem 1.2 of 
5l| lor see 


3. It can be proven directly by a slight modihcation to the proofs 


in 


5^), in which case it comes ont natnrally as a mnltidimensional vari¬ 


ant, where the exponent '^{E) is connected with decay of generalized eigenfnnctions 
in the mnltidimensional space. We include this theorem here, showing that it follows 
from Theorem |^, mainly to illustrate the connection of Theorem with this type 
of results. 


It is often convenient to filter scaling behaviors in terms of explicitly defined scal¬ 
ing exponents and to formulate relations between scaling behaviors as inequalities 
between such exponents. It is thus natural in our context to define 


(3 = lim limsup 

5^0 T^oo 


log (inf [L I (||e > ^}) 

logT 


B = lim lim inf 

— 5^0 T^oo 


log (inf [L I (i|e > <5}) 

logT 


(3 and [3 are the upper and lower spreading rates associated with the slowest spreading 
portion of the wavepacket. Local exponents for the asymptotic scaling behavior of 
the scales L^{E) are given by 


A(E) = limsup 

loge-i 


ME) = lim inf 


log4*(-E) 

log 6“^ 


Theorem 1.1 immediately implies 


(3 < fi'^-essinf X{E). 

Moreover, we see that in problems where the solutions uo{n, E) behave nicely enough 
to ensure that ||mo(77, E)\\L^/\\uo{n, E)\\l^ is small whenever L 1 /L 2 is small. Theorem 


1.3 would imply 


inf A(i?) < (3. 

These inequalities are particularly interesting in cases where the problem exhibits 
nicely scaling power law behaviors so that uo{n,E) behaves as described above and 
also //■'■-ess inf A(i?) = //■'■-ess inf A(i?). If this happens, we see that we also have 
f3 = f3 and the spreading rate of the slowest spreading portion of the wavepacket is 
precisely given by the slowest growth rate of the scales L^E) with respect to the 
spectral measure //■'■. 

Our second main goal in this paper is to apply the new dynamical upper bounds 
to the Fibonacci Hamiltonian, the most studied of all one-dimensional models of 
quasicrystals. For this application we need a whole-line (namely, version of 

Theorem O- To formulate this whole-line version, note that the scales [E) and 
L~{E) can be defined in a way directly analogous to (|1.5|), involving the same 

kind of sums but taken over the negative half-line. Notice that the negative half-axis 
in our setting is (...,—1,0), so that the summation in the analog of definition (p..5|) 
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will now start (or rather end) at 0 instead of —1. Let /i be the spectral measure of 
the whole-line operator H corresponding to (5i, dehned by 


M(z) ^ ((H - 

Jr X- z 


Theorem 1.5. Let H be a self-adjoint operator of the form (0) on Let the 

characteristic scales 
we have 


L^_i(-E) he defined by ( |1.5| ). 


Then for any T > 0, Li, L 2 > 1, 


(1.13) (||e-*^*5i||i^,i,)r > Cfi ({E I L-_,{E) < L, and L+_,(E) < L 2 }) , 

where C is a universal positive constant. If Li, L 2 > 2, the same result holds if 
L^_i{E) is replaced in the statement with L'^_i{E) defined by (O). 


The Fibonacci Hamiltonian Hx is dehned by 

{Hxu){n) = u{n -l- 1) -l- u(n — 1) -l- AH (n)u(n), 

where V(n) = [(n-(- l)a;J — [na;J and co = (V^— l)/2 is the golden mean. It is known 
P, ^ that for every A > 0, ffx has purely singular continuous spectrum, and 
its spectrum (as a set) is a Cantor set of zero Lebesgue measure. Lower 

0 


moreover 

bounds on wavepacket spreading rates for ffx were recently shown in 


We are 


going to show both upper and lower bounds for the dynamics of Hx which imply that 
the spreading rate is intermediate between ballistic (~ T at time T) and localized 
(~ T°). To the best of our knowledge, this is the hrst proof of such behavior in an 
explicit model of this type. 


Theorem 1.6. Let Hx be the Fibonacci Hamiltonian. Then 

i) There exists a constant G > 0 such that for sufficiently large A, 

(1-14) ( 11 ^ ^ VT > 0, 

where pi(A) = C'i(log A)“^(l -|- 0((A log A)“^)). 

ii) Given any g > 0, there exists a constant Gg such that for every coupling A > 0 
and every time T > 0, 

(1-15) ^||e ^ ^illcgTP2(^),CgTP2(^)^.^ — 

where the positive function p 2 {X) satisfies p 2 {X) = C 2 (log A)“^(l-|-0((log A)“^)) 
for large X. 


Remarks. 1. The hrst part of Theorem |1.6| is an upper bound on wavepacket spread¬ 
ing for large coupling, which shows that on the average, there is a portion of the 
wavepacket remaining in a ball of the radius ypi(A) T. The second part of 

the Theorem provides a lower bound on wavepacket spreading, showing that on the 
average, only an arbitrarily small part of the wavepacket remains in a ball of radius 
~ at time T. The constants pi(A) and P 2 {X) have the same (logarithmic in A) 

asymptotic behavior for large coupling, up to a constant factor in front of the main 
term. This implies that the logarithmic law is precise and cannot be improved in the 
estimates (|1.1^ ), ( |1.15|) . 
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2. It will follow from the proof below that A > 8 is sufficient to get a nontrivial 
upper bound (pi(A) < 1 in (|1.14|) ). This range is not optimal and can be improved 
by additional technical effort; however our current methods do not allow to extend 
the bound all the way to A = 0. Moreover, if A is so small that pi(A) > 1, then there 
is little interest since this would constitute a ballistic (or worse) bound. 

3. We remark that Sinai, in a recent paper [^, studied anomalous transport (in 
terms of moments of the position operator) for an almost periodic potential in a dif¬ 
ferent setting (in fact, in a pure point regime). 


Although we only treat discrete operators in this paper, results parallel to Theo¬ 
rems o, PI, 0,0,0 also hold in continuous settings by direct translation of 
the arguments given here. 

The rest of this paper is organized as follows. In Section 2 we prove some auxiliary 
facts relating the scale L'^{E) and the m function. In Section 3 we derive the upper 
bounds on wavepacket spreading. In particular, the proofs of Theorem |p. Theo¬ 
rem IP and Corollary |1.2| appear there. In Section 4 we show the new lower bounds. 
Theorems |1.3| and 0. In Sections 5 and 6 we treat the Fibonacci Hamiltonian and 
prove Theorem |1.6| . Finally, In Section 7 we discuss how our results are applicable to 
multidimensional problems. 


2. Bounds on the m function 


We are going to prove a series of auxiliary estimates relating the behavior of 
solutions Mo, 7 r /2 and the half-line m functions. The line of the argument follows 
Jitomirskaya-Last extension of subordinacy theory |p|, but there will be an essential 
modihcation that will be crucial for the derivation of the dynamical criteria. 

Recall the following simple example of variation of parameters which will prove 
extremely useful: 


Lemma 2.1. Suppose w{n, z),w{n, E) solve o for spectral parameters z, E re¬ 
spectively. Assume that tc(0, z) = tc(0, E) and tc(l, z) = tc(l, E). Then for n >0, 


n 

(2.1) t(;(n, z) = w{n, E) -\- {z — E) A;(n, m; E)w{m, z) 

m=l 


where 

k{n, m; E) = uo{n, E)uT,/ 2 {m, E) - UT,i 2 {n, E)uo{m, E) 

and with the convention that the sum is zero if the upper limit is less than the lower 
one. 


Proof. It is a direct computation to check that the right-hand side satisfies (|1.3|) with 
the spectral parameter z. Hence, on both sides we have solutions of (|1.3|) which 
coincide at the sites 0,1. This implies the equality. □ 


Equation (|2.1| ) has the form 

w{n, z) = w{n, E) + {z — E){K{E)w){n, z), 
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where K{E) is an integral operator with the kernel k{n,m] E). As introdnced in 
( |1.4| ), II ■ IIi dehnes a norm on an [LJ + 1-dimensional Hilbert space (except when L 
is an integer, in which case the dimension is L). The corresponding inner prodnct is 


[L\ 


(2.2) (01, 02)L — 4>i{n)4>2{n) + {L — [Lj)0i(LLj -|- 1)02( [LJ -|- 1). 


n=l 


We wish to estimate the norm of the operator K{E) acting in this space: 


Lemma 2.2. The Hilbert-Schmidt norm of K{E) is given by 

(2.3) ll|i^(i?)llli=ll«o||ill«./ 2 lli-|(«o,«./ 2 )nP 

(2.4) = msix\\ue{n,E)\\l ■ min \\ug{n, E)\\l. 

9 9 

where are the solutions at energy E. In particular, this gives an upper bound 

on the operator norm of K. 


Proof. Recall that the Hilbert-Schmidt norm of an integral operator is eqnal to the 
norm of its kernel. For integer L, 

L 

|||A'(B)|||I= ^ |fc(n,mT)r = ll«ollill«./2fi-|{«o,«,/2)Lr 

n,m=l 

can be shown fairly easily. Non-integer L merely reqnire a more lengthy compntation. 
To show the second eqnality, consider the 2x2 matrix 

Q^ \ hoWl {Uo,U^/2)l 

('^7r/25 ^o)-L ||^7r/2||7;, 

It is self-adjoint and positive, and if eg = (cos 6 *, sin^), then 

egQef = ||m 6 i||l- 

The right-hand side of (^1^) is the determinant of Q and so the prodnct of the 
eigenvalnes of Q. By the minimax characterization, these eigenvalnes are given by 
the factors in (|2.4| ). □ 


Remark. Notice that if the maximnm of ||M 6 i(n, i?) ||l is achieved at 6 *o, the minimnm 
is achieved at the orthogonal bonndary condition 6^ (n^x satisfies M 0 x(O, E) = cos 6 ^ 0 , 
Ug^{l,E) = — sin 6 *o). This is becanse the eigenvectors of the self-adjoint matrix Q 
are orthogonal. 

From (P^) we know that |||iF(F^)|||L is strictly increasing. Hence e|||iF(i?)|||L = 1 
determines L as a fnnction of e and E. As in the introdnction, we denote this length 
scale by Lf{E). The main resnlt of this section is an estimate for the norm of the 
Weyl solntion over this length scale. Before stating and proving this result, let us 
recall some facts about the Weyl theory of Jacobi matrices (see, e.g., [^): For each 
2 ; G C\M, there exists a unique solution of (|1.3|) which is square summable and obeys 
M+(0, z) = 1. One may write this solution in the form 


n+(n, z) = M^/ 2 (n, z) - m+{z)uo{n, z), 
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and so define the Weyl m-fnnction m^{z). As we will discuss a little more in the next 
section, m^{z) captures all of the spectral information about For the moment, 
however, we merely need the observation that 

CO 

(2.5) a(0)lmm_|_(z) = Imz E \u+{n,z)\‘^, 

n=l 

which follows from summation by parts. Recall that we assume a(0) = 1. 

Theorem 2.3. Fix G M and e > 0, then we have 

( 2 . 6 ) 


2-v/3<ll“°llyl""P + ")l<2 + V3. 


Moreover, 

(2.7) 


flk+IlP > 


Il'^7r/2||L+ 
\m+{E + ie)\ 


4e||Mo|lL+||«^/2|lL+ 

In both formulae, uo,Unf 2 are the solutions at energy E 


> ^-j^lmm+{E + ie). 


Remark. The bounds (|2.6|) are very similar to those proved by Jitomirskaya-Last p2 
The difference is that the scale Lf is given in by the condition 

^II'^o|Il+I l'*l7r/2||L+ = 1- 
The scale Lf which we dehne by 

(2.8) |||A:|||^+=max||ne(n,E)||2+ ■ m:m\\ue{n,E)\\l+ = e-^ 


e e 

is larger (compare (p.3|) ), and the constants appearing in 
^| . The scale defined by 


m 


are sharper than those 
might be less convenient for dimensional spectral 
analysis, since its dehnition is more involved. However, an improvement contained 
in (|2.8|) will be quite crucial for the proof of dynamical criteria in the next section. 

Proof. To shorten the formulae which will follow, we introduce 

{uq, UT^jf) 1^+ 




and C = \/l — /3^. 


ll'*^o||L+ll'^7r/2|lL+ 

Notice that Cauchy-Schwarz says (^ > 0. At times we shall also write for m^{E + 

ie). Notice that by Lemma (pT3|) and (pTS]), the scale Lf{E) given by ([T7^) is 
chosen exactly in a way to ensure that ||iF(i?)||L < e~^ if L < Lf{E). Lemmas |2.1|JT2 
and the definition of Lf combine to show that 

(2.9) 4\\u+{n, E + ie)\\l+ > ||M^/ 2 (n,^) - m+{E + ie)uo{n, E)\\l+. 

Using ( p.5|) , we see that 

4e“^Imm+ > ^\u+\W+ > \\u.,,/ 2 \\\+ + |m+p ■ ||mo||^+ - 2Rem+(Mo,M7r/2)L+- 
The definitions of (3X and Lf are such that eCII'Wo|lL+Il'*i 7 r/ 2 ||L+ = 1- Thus, 


(2.10) Imm+ > e||M+||^+ > 


\u. 


r/2||j 


4C I ll%lh 


+ \m+{E + ie)\ 


2 II^oIIl+ 

Il'^7r/2||L+ 


— 2/3Rem+ 
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This equation implies both claims of the theorem. To prove ( p.6|) , notice that | ± /? + 
zCl = 1 and so 

—2|/3Re (m+)| — 4CImm_,_ > —2|m+| — 2Re [(±/9 + iOm^] > —4|m+|, 
where the sign is chosen to make ±/3Re (m+) > 0. Applying this to ( p.lOD , leads to 


II«oIIl+ 


\m^{E + A)p — 4:\m^{E + ie)\ + 


I'^7r/2||L+ 


Il'^7r/2||L+ ir'^0|lL+ 

which is equivalent to ( f^.bl) . To prove ([^.7|), use CT) to obtain 


< 0 , 


e\\u+{n, E + ie)\\l+ > 


\m+\ ( ||M7r/2|lL+ , ll“ollL+l"i+l 


4C l|m+|i|no|L 


+ 


l’^7r/2|| j 


2/3- 


m_i 


> 


l"^+lio OI/PII C\m+\ 


4C 


-(2-2|/3|) = 


2(1 + |/3|) 


> \C\m+{E + ie)\. 


(We used x + x ^ > 2 in the second step.) This proves the left-hand inequality in 
because eCII'^^o|lL+||'W 7 r/ 2 ilL+ ~ complete the proof, we need to show the 

right-hand inequality. From ( ^.9| ) and ( |2.5|) we infer that 

4Imm+ > e(Imm+)^||Mo||^+, 


and from (p.6|) , 

\\uo\\l+ > {2 - V3)\m+{E + ie)\-^\\uo\\L+\\u^/ 2 \\L+- 

Combining these two gives the right-hand side of (|2.7|). 


□ 


The following variant of the inequality ( |2.7]) in Theorem |2]^ will prove useful in 
the next section. 


Proposition 2.4. Suppose E, E' G 
universal constant C, so that 


e > 0 and \E — E'\ < e. Then there exists a 


( 2 . 11 ) 


\u+{n, E' -t- - ClTam+{E' + ie). 


Remark. It is important that while the bound in (p.ll|) involves m+ and the solution 
M+ at the energy E' + ie, the scale Lf is defined at the energy E. 


Proof. This result is a direct corollary of the proof of Theorem |2.3| . All statements of 
Theorem ^]^hold with E + ie replaced by E' + ie and adjusted constants. Specifically, 
the only change introduced by replacing M+(n, E + ie) by n+(n, E' + ie), \E — E'\ < e, 
is a change in the constant in (|2.9| ) from 4 to (1 -|- a/ 2)^ < 6. Following through the 
proof with the constant 6 shows that C may be chosen to be (3 — 2v^)/36. □ 



















BOUNDS ON WAVEPACKET SPREADING 


13 


3. An upper bound on wavepacket spreading 

We begin by relating the dynamical quantity we need to estimate to a solution of 
equation (|1.3|). As before, 

u+(n, z) = m^/2(r-, z ) - m+(2;)Mo(n, z) 

denotes the unique solution for z G C \ M with m+(0, z ) = 1. 

Lemma 3.1. For every z G C \ M and n > 1, 

(3.1) {{H+- z)~^6i,6n) =-u+{n,z). 

Proof. It’s easy to see from the dehnition of the resolvent of that {{H^—z)~^5i, dn), 
n > 1, is an solution of ( |1.3| ). It thus follow from the uniqueness of the Weyl 

solution that it must be a multiple of u+(n,z). The Lemma is thus implied by the 
dehnition of m+(z) and our normalization of u+(n, z). □ 


Lemma 3.2. For any T > 0, 

(3,2) \u^(n,E + f)fdE. 

K 


Proof. Recall that 




The right-hand side in the above equality is a constant times the square of the norm 
of the function f where the complex measure /ii^„ is dehned by 

Jr 

The Fourier transform of this function is equal to 


X — z 


j_ j dizUE') _ _ B s,.) = J-uAn. E + -] 

271J {E-E') + {i/T) 27r'^ + ^ 2?! 


by (0. Now (|3.2|) follows from Parseval’s equality for the Fourier transform. □ 
Lemma |3.2| implies, 

(3.3) ^ F j Wu^n, E + ^)\\ldE. 

R 


Given a set S' C M, let us denote by S'^ the e-neighborhood of the set S. We now 
complete the proof of the half-line dynamical bound. 

Proof of Theorem O- Let 

S = {E\L+_,{E) < L}. 




14 


ROWAN KILLIP, ALEXANDER KISELEV, AND YORAM LAST 


By Lemma to prove Theorem 0 it suffices to bound the right-hand side of (p.3|) 
from below by Hence, it suffices to estimate from below the integral 


e J \\u+{n,E + ie)\\ldE, 
s. 


where L > L'^{E) for E & S (we think of e = T ^). By the definition of for every 
E' G S'e there is G S' such that \E — E'\ < e. Applying Proposition we obtain 


(3.4) 


||M+(n, E ie) 


\idE>C 


lmm_^{E -|- ie) dE. 


By the Fubini theorem, the left-hand side in (|3.4|) is equal to 


C / dfij^{x) 


edE 


{x - E)2 + e2 


> (7 / dn+{x) 


edE’ 


(E')2 + e2 


> Cf^+{S). 


□ 


The following Lemma shows that Theorem 0 implies Corollary |1.2| . 

Lemma 3.3. For every G M and L >2, 

(3.5) 4||4(l,i5)-‘f (max||t.„(n,B)||l) (mm ||«»(n, B)||i) > ||4(i5)Hi. 

Proof. A direct computation using the dehnition of <F(n, E) shows that for any 9, 
\\^{n,E)\\‘^ < {u 0 {n,E))^ + {ue{n + l,E))'^ + {u 0 ±{n,E))‘^ + {ug±{n + l,E))‘^. 


Thus, 

\\^{E)\\l < 2{\\ue{n,E)\\l + \\ue±{n,E)\\l) < 4 (^max ||ne(n, E)||i) . 


For L > 2, we have for any 6 , ||m6»IU — ^ 11^(1)-^) 

(|3^ ) holds. 


-iii-i 


and we thus see that 
□ 


Proof of Corollary \1.2[ . Since Lf{E) was defined in (|1.7|) by 

we see from o and Lemma that Lf{E) > Lf{E) whenever Lf{E) > 2. 
Therefore, Corollary |1.2| follows immediately from Theorem |1.1| . □ 


It remains to consider the whole line case. Denote by M_(n, z), z G C\M, a solution 
of ( p..3|) which belongs to = {0, —1,...,}, and is normalized by 


M_(n, z) = u.^i 2 {n, z) -h m-(z)uo(n, z). 
Recall the notation H for the operator defined by (|1.2|) on C{Z). 
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Proof of Theorem \1.5^ . As in Lemma |3]^, we can show 

(lk-“%llLi,jT = ^ j \\a(l,n,E+'- 

'J M. 

where 

G(l, n,z) = (fH - 6n). 

One can easily verify that 

G(l,n,z) = | n>l 


-^- It { Ti 77 ^ 1 

a(0)(m+{2)+m_{2)) “'-OU ^ 


In particular, 


G{l,l,z) = M{z) = 


-m+{z)m{z) 


a{0){m+{z) + m_{z)) ’ 
and the spectral measure fi of H corresponding to the vector (5i satishes 

djji^x) 


M{z) = 


X — z 


Now define 


S = {E\ L+_,{E) < La and Lf_,{E) < Lj . 


Repeating the same arguments that led us to the proof of Theorem o, we get 
(6 = T-1) 


e\\G{l,n,E + te)\W^^>C 

i 

and similarly 

e\\G{l,n,E + te)\\l^,>G i 


\m^{E + ie)\" 


\m+ {E + ie) + m_ {E + ie) \ 


-Im m+{E + ie) dE 


\m+{E + ie)\" 


|777,+ (E + ie) + m-{E + 7e)p 


lmm-{E + ie) dE. 


Combining these two inequalities, we get 

e\\G{l,n,E + ie)\\l^^L 2 >C [ lmM{E + ie)dE > C/i(L), 


similarly to the last step in the proof of Theorem 0. 


□ 


4. A LOWER BOUND ON WAVEPACKET SPREADING 


Before treating our main application, the Fibonacci Hamiltonian, we pause to prove 
the general lower bound on dynamics given by Theorem 0. 
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Lemma 4.1. Let Ps be a spectral projection of the operator on a measurable set 
S. Then 


(4.1) 


(e-*^+‘Ps5i) in) = / e-^^^Uo{n,E)dfi+{E). 


□ 


Proof. Since is a cyclic vector for it follows from the spectral theorem that 
is unitarily equivalent to multiplication by the parameter on d/i’*'). UQ{n,E), 

n > 1, is known (see, e.g., P) to be the representation of 6n in this space. 

Proof of Theorem E3 Notice that according to Lemma E3 

oo 

2 


< 


T 


e / e^^^UQ{n,E) dpP{E) / e ^^'^Uo{n, E') dp,'^{E') dt 


n<L g 


r r 4T“^ 

^ / / uo{n, E')uo{n, E) 


n<L g g 


{E - E'Y + 4T-2 


— d/l+(p)d/l+(p') 


(4.2) 


< 


4T 


-2 


-dp+(P') I ||«o(n,P)||id/i+(P). 


{E - E'Y + 4T-2 

s \s 

We used Cauchy-Schwartz in the last step with respect to the product measure 
dpP{E)dpP{E'). Notice that the term in the brackets in (^^) does not exceed (setting 
e = T-i) 

2eImm+(P + 2ie) < 4eImm+(P + ie). 

By the estimates ( ^(^ and (|2.7|) of Theorem |^, we have 
(4.3) 

I JC-i -MI ( I7MI2 ^20, /2;,fI“™+(.E + !e)||tlo|lt+||M,/2lllJ , 

elmm+(h; + ie)\\uo[n, E)\\^+ < (2 + V3)-^-tv;-ttt-^ < C. 


\m+{E + ie)\ 
ho(n,P)||? 




l|wo(ri,p)ih 


Combining (14.2|) and (O), we obtain 

{\\e-^^+^Ps54l)T<C j 

s 

which is exactly what we wanted to show. □ 

Proof of Theorem O- By the extension of subordinacy theory due to Jitomirskaya 
and Last 


lim inf 


holli-“ 


= Oy»P“/i+(P) = cx). 


L^oo ||'Uvr/2||j 

Therefore, for a.e. E, we have with some C'i(P) > 0 


2-a{E) 
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By (^), the definition ( |1.5| ) of L+ and relation (|^) , 


-1 


— II^o||l+I| ll7r/2||i+ Cl (£^)||rto 


2 

I 


Therefore, 

(4.6) ||'t.„||2+ > 

On the other hand, by the definition of '^{E), we have 

Iholli < C2(i!)L’W, 


Denote 


5 = {E\'q{E) > h} 

(recall rj{E) = a{E)/'^{E)). Let S'! C S' be the set snch that for E G Si, Ci{E)C 2 {E) < 
C with some uniform constant C. Clearly, Si may be chosen so that ||T5\5^(5i||^ = 

we find (with C denoting different universal constants in different places) 

(4.6) (l|e-“+Vs,6,||i)T <C [ 

J \\uo{n,EW + 

Si 

Si 

Choosing L = CgT^, we can make the left hand side in ([4.61) arbitrarily small by 
changing Cg. It remains to observe that 

(l|e-*^+*5i|li)T < ((||e-*^+*P5i5ilU + ||e-*^+*(/ - Psi)5i||l)^)t 

We claim that by choosing Si and then Cg the last expression can be made smaller 
than 1 — fi'^{S) + g for any g > 0. Indeed, first choose Si so that 

\\{I-Ps,)6ir<l-f.HS) + {g/2). 



IJ.~^{S \ Si) is as small as we want by adjusting C. Then from (p..9|) of Theorem 1.3 


Then choose Cg so that 


|e-*^+*P5i5i|li)T<(l?/5) 


(assuming g is small). This completes the proof, demonstrating that the part of the 
wavepacket corresponding to the energies in S leaves (on the average) a ball of the 
size ~ at time T. □ 
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5. The Fibonacci Hamiltonian: bounds on traces 


Let us recall from the Introduction that the Fibonacci Hamiltonian is the discrete 
Schrodinger operator 

(5.1) [Hxu] (n) = u{n + 1) + u(n — 1) + XV (n)u{n) 

acting on with the Fibonacci potential V{n) = \ {n + l)a;J — \nuj\. Here uj is 

the golden ratio, (\/5 — l)/2. The most important property of this potential is the 
substitution rule p8| , 

(5.2) V{qk + n) = V{n) for n = 1,..., and fc > 3 


where qk denote the Fibonacci numbers, go = 1) <?i = 1) <?fc 
For /c > 1 we dehne by <hfc the transfer matrix <h(gfc, E) : 


u{qk + 1 ) 


«(!)■ 

1 - 

— 

“( 0 ). 


Qk-l + 


for all u{n) satisfying H\u = Eu. It is convenient to make the additional dehnitions 


$-1 = 


1 -A 
0 1 


Tn — 


E -1 
1 0 


From the substitution rule it follows that d>fc+i = <hfc_id>fc, k > 0 [21]. Let us denote 


Xk{E) = Tr<hfc. With more work, one obtains the trace map and trace invariant [^|: 


(5.3) Xk+i = XkXk-i - Xk-2 

(5.4) Xfc+i +xl + xl_^ - Xk+iXkXk-i = 4 + 


We denote by ak the spectrum of the periodic potential with period qk coinciding with 
V(n) for n = 1,..., qk- By the Bloch theorem, ak is a set of intervals (bands) for which 
Xk{E) G [—2,2] (see, e.g. [@). Moreover, in each band Xk{E) varies monotonically 
in [—2, 2] and takes values ±2 at the ends. Due to the relations (|5.3|) , (|5.4|) , the traces 
are among the most convenient objects of study in the Fibonacci model. Since Xk 
describe the spectrum of periodic approximants, there is also a natural relation to 
the spectrum of the limiting Fibonacci Hamiltonian. The following properties of the 
traces and their relation to the spectrum of H\ are quite useful: 


Proposition 5.1. i) The spectrum of H\ coincides with the set of energies for 
which the sequence Xk{E) is bounded. 

ii) If \xk{E)\ > 2, Xk+i{E) > 2 for some fc, then the sequence Xn{E) is unbounded. 
hi) If A > 4, there cannot exist E,k such that \xk{E)\ < 2, \xk+i{E)\ < 2 and 
\xk+ 2 {E)\ < 2. 


Proof. The hrst and second statements have been proved by Suto 
statement is a direct consequence of the trace invariant (|5.4|). 


38 


The third 

□ 


Our goal in this section is to prove a lower bound on the derivative [xj.(£')[ for 
energies in the spectrum of H\. In the next section we show that this bound translates 
directly into the lower bound on the growth of \\^{n,E)\\f that we need for the 
upper bounds on dynamics. The results of this section are contained, in a somewhat 
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k + 2 
k + 1 
k 


k-1 


Figure 1. Types of bands in cr^ : left, a type A band; right, a type B band. 


different form, in the preprint of Raymond |^, who nsed them to derive an upper 


bound on the Hausdorff dimension of the spectrum (which alone does not imply any 
upper bounds on dynamics, see ||^). For the sake of completeness, we present here 


a simplihed version of the argument given in |^. The main result of this section is 


Proposition 5.2. Assume that the coupling A is sufficiently large (\ > 8 will do). 
Then for every E in the spectrum of H\, the derivative of the trace Xk{E) satisfies 

(6.6) |i;(B)| > 

where .^(A) > 1 and 

e(A) = A(l + 0(A-i)) 

in the large coupling regime. 


Remark. We do not attempt to get the optimal range of values of A. Instead, we opt 
for the clarity of exposition and freely assume that A is large enough. The arguments 
presented below lead to nontrivial dynamical bounds for A > 8. With more technical 
effort, this value can be reduced, but remains far from zero. 

Definition. We call a band Ja, C cjfc a type A band if R C ak-i (and consequently 
by Proposition 04 n (cTfc+i U ak- 2 ) = 0 )- We call a band Ik C Ck a. type B band if 
h C ak -2 (and so 4 n cjfc-i = 0 ). 

See Figure 1 for an illustration. Notice that by the dehnition of <h_i and $ 0 , 
cr_i = M, (To = [—2, 2] and ai = [A — 2, A + 2]. Hence, (Tq consists of one band of type 
A, and ai consists of one band of type B. The following Lemma allows to determine 
inductively the structure of the spectrum of II\. 

Lemma 5.3. Assume that A > 4. Then for any k > 0 : 

i) Every type A band Ik C ak contains exactly one type B band I k +2 C a k+ 2 , and 
no other bands from ak+i, ak+ 2 - 

ii) Every type B band Ik C ak contains exactly one type A band R+i C ak+i and 
two type B bands from ak+ 2 , positioned around Ik+i- 

Proof. Consider a type A band Ik (Z ak- By dehnition, Ik C ak-i, and so |xa;+i| > 2 on 
Ik by Proposition b.l|iii). This implies IkHak+i = 0 . On 4 , Xk changes monotonically 
from —2 to 2, in particular, there is a unique Ei G h such that Xk{Ei) = 0. By the 
trace map (|5.8|), we have \xk+ 2 {Ei)\ = |a;A;-i(F'i)| < 2, so 4 0 ak +2 7 ^ 0 - Notice also 
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that when Xk = ± 2 , \xk+ 2 \ > 2 |a;A;+i| — |a;A;_i| > 2 , so all possible bands of ak +2 
intersecting lie inside Ik- Moreover, in each band Ik +2 C Ik of crfc+ 2 , Xk +2 changes 
from —2 to 2 , and so by intermediate value theorem there exists E 2 G Ik +2 such that 
Xk+2{E2) = -Xk-i{E2). Then by (|^) , Xk{E2)xk+i{E2) = 0, and since \xk+i{E)\ > 2 
in Ik, Xk{E 2 ) = 0. Hence, every band Ik +2 C h contains an energy where Xk = 0; by 
monotonicity, there is only one such band. 

Now consider a type B band Ik C (Tk- 2 - When Xk{E) = 0, we have |xfc+i| = \xk- 2 \ < 
2 , and so Ik H ak+i 7 ^ 0- As in the above argument, |a:fc+i| > 2 when \xk\ = 2 , so 
all bands of Uk+i intersecting Ik lie strictly inside. Moreover, by same argument as 
above, any band Ik+i C ak+i inside Ik must contain an energy E where Xk{E) = 0. 
Hence, there is a unique such band Ik+i C Ik- Next consider ak+ 2 - Iterating the trace 
map, we hnd Xk +2 = — l)xk-i — Xk- 2 Xk- When Xk = ±1, \xk+ 2 \ < 2. Also, if 

Xk = ±2, \xk+ 2 \ > 3|a:fc_i| — 2\xk-2\ > 2. Therefore, there are at least two bands of 
ak +2 lying strictly inside Ik to the right and to the left from Ik+i (bands of ak +2 C h 
cannot intersect Ik+i by iii) of Proposition |5.1|) . It remains to show that there are 
only two such bands. Using (|5.3|) , it is easy to verify that 

(5.6) {xk ± l){xk+2 ± Xk- 2 ) = {xl - l)(a;A:+i ± Xk-i)- 

Consider a band Ik +2 C ak +2 lying in Ik- For E G Ik+ 2 ^ Xk+i{E), Xk-i{E) have hxed 
signs (both satisfy \xk±i{E)\ > 2 there). Pick a sign in ( |5.6| ) so that Xk+i ±Xk-i 7 ^ 0. 
By the intermediate value theorem, there exists the energy E G Ik +2 where Xk +2 
Xk -2 = 0. At this energy, we must have x^ — 1 = 0- But there are only two energies 
in Ik where Xk{E) = ±1. Hence, there are at most two bands of ak +2 in Jfc. □ 

For the proof of the key lemma, we need the following auxiliary result. 


Lemma 5.4. Let the functions f±{x,y,X) be defined by 
(5.7) f+{x,y,X) = I xy ± 

Assume that A > 4. Then for —2 < x,y < 2, 


( 6 . 8 ) 


dU 


dx 


A,y,X) 


df+ 


dy 


x,y,X) 


< 1 . 


Proof Since f_{x,y,X) = -f+{x,-y,X), and f+{x,y,X) = U{y,x,X), it suffices to 
show the bound for 


9f+, .. 1 I 


x{4 - y^) 


^4A2 + (4-a;2)(4-j/2) 


Notice that 


d^f+ 1/ SxyX'^+ xy{A-x^){A-y^y 

(x, I/, A) = - 1 + - ttwttw 1 > 0 


dxdy ^2 
for A > 4 and |a;|, \y\ < 2. Thus, 


max 

hlUI<2 


dx 


{x,y,X) 


(4A2 + (4 — a;2)(4 — i/^))^/^ 
df. 


= max 

\x\<2,y=±2 


dx 


{x,y,X) 


= 1 . 


□ 
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The Proposition |5.2| follows immediately from Lemma and 


Lemma 5.5. Assume that the coupling A is sufficiently large. Then there exists a 
number > 1 such that, given any (type-A) band Ik+i C ak+i lying in the band 
Ik (T ak, we have 

4+1 (D 

4(D 



for E G Ik+i- Similarly, for any (type-B) band Ik +2 C <7^+2 lying in Ik, 


x'k+2i.E) 


>e(A) 


for E G Ik+ 2 - Moreover, = A(1 + 0(A ^)) in the large coupling regime. 


Proof. The proof is by induction. The induction is well-founded since x'_i = 0, x'q = 1, 
x'i = 1. There are three cases to consider. 

1. First, consider the case where Ik+i C Uk+i lies in Ik C ak- Then by Proposition |5.1| 
Ik n = 0, and Ik C crk- 2 - Differentiating the equality Xk+i = XkXk-i — Xk -2 and 
dividing by we hnd 


(5.9) 


x't+i _ ^ 

— Xk-1 + 


X 


k-2 


XV 


XV 


XV 


From the trace invariant (El 

2 , 2,2 + I \2 

Xk + X^_-^ P X^_2 - XkXk-lXk-2 = 4: + X 


we get 

Xk-i = f±{xk, Xk-2, A) = i (^XkXk-2 ± y^4A2 -i- (4 - a;^_2)(4 - xD'J > A - 2 

since \xk\, \xk- 2 \ < 2. Therefore this term in ( p.9|) is large. To estimate the remaining 
terms, notice that 

Xk—l (,Xk, Xk—2-i ^'jXk T (,Xk, Xk—2-i ^')Xk—2 

(where ± means that either plus or minus may occur). Applying Lemma |5.4| , we 
derive from ( p.9|) that 

P, . . P, „ 

- \Xk\- 


^'k+1 

> A-2-3 

x'k-2 

^k 


^‘k 


By the induction assumption, \x'j^_ 2 /x'k\ < 'C(A) so the induction step holds true 
provided that 

e(A)<A-4-3e(A)-\ 


The maximal ^(A) we can take under this condition is 

5(A) = 1 ((A - 4) + V(A-4)2-12) = A + 0(1), 
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Next, consider Ik +2 C < 7^+2 lying in 1^. Here we have to distingnish two scenarios. 
2. 4+2 n ak-i = 0, 4+2 C crfc- 2 . Given 


/ f , f / 

- ^y» ^y» I ryt ry* ry> 

'^k -\-2 '^k+l'^k ' '^k+l'^k 15 

/ / I / / 

/7> - ry> /y> \ ry> ry> ry> 

■^k+1 •^k‘^k—1 W •I'k'kk—l ■^k-2 


we find 


/^' ry»' ry>' 

'^k+2 o I /' 2 ‘‘'fc-2 

-^ = 2Xk+l - Xk-2 + (Tfc - 1)—;-Tfc—. 

/V>' V 11/ / „/ „/ 

•^k -^k -^k 

Similarly to the previons argnment, 

\xk+i\ = \ f±{xk+ 2 ,Xk, A)| > A - 2. 

Also, 

Xk—X \Xk) Xk—2i ^jXk \Xki Xk—2i ^)Xk—2' 


Using Lemma |^ , the bonnds \xk\, \xk- 2 \ < 2, and the indnction assnmption, we 
arrive at 


x‘ 


k+2 


X', 


>2A-9-5^(A) 


-1 


Solving the qnadratic ineqnality, we find that the maximal ^(A) for which indnction 
step goes throngh in this case is 

i{X) = i (^(2A - 9) + V(2A - 9)2 - 20 ) = 2A + 0(1). 

3. 4+2 C crfc_i, 4+2 n (Tk -2 = 0- In this case we obtain 


X' 


k-\-2 I ft-i 

— = Xk+l + Xk -7 


^fc+1 ^k-1 


X 


^k ^k 


As before, x^+i > A — 2. Also 

df± 


9f± 

^fc+l {Xk^ Xk—\i A)x^ A {Xk ) Xk—\i X^Xk—i- 


This leads to 


X 


k+2 


X', 


> A-4-3^(A) 


-1 


Hence the indnction step works with the same ^(A) as in the first case. 

From Lemma |5.3| it follows that we have considered all possible sitnations. 


□ 


Remark. In particnlar, the given argnment shows that Lemma b.5| (and hence Propo¬ 
sition |5.2|) holds with ^(A) > 1 for A > 8. 
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6. The Fibonacci Hamiltonian: dynamical bounds 


Our first goal is to relate the derivatives of traces with the information required by 
the dynamical criterion of Theorem O- The following formula is well-known (see, 

e.g., 0) 

qk 

^kiE) = {^u^/2{qk, E)ul{n, E) - + 1, E)ul/^{n, E) + 

n=l 

+ [u^/2{qk + 1,F) - Uo{qk,E)\uo{n,E)u^i2{n,E)^ 

and can be derived using variation of parameters. Therefore, 

<(B)<4||4(n,i5)||=,_H 

(one can show that x'^{E) < 4(A -|- 0(l))||<h(n, F)||q^_,_^ for E in the spectrum of Hx, 
but we choose to avoid the technicalities). Hence, Proposition |5.2| implies 

(6.1) ii4(>i.B)ir ' ^ 

As oj is the golden ratio {y/E — l)/2, 

qk = [uj~^ - {-ujY]/y/E 

for /c > 1 (see, e.g., [^]). So for large k, qk ~ uj~^lx/E. Consider qk "E L "E qk+i- 
Since qk+i/qk < 2 for all /c, (|6.1|) implies 


L+. > 


( 6 . 2 ) 


||4>(n,B)iri > > C'L<‘ witli Cl = 


Later we will derive an upper bound on this quantity of the form CL^'^ with C ,2 also 
depending logarithmically on A. 

Now we complete the proof of the upper dynamical bound in Theorem |1.6|. 


Proof. By the symmetry of Fibonacci potential (F(—n) = F(?7, — 1) for n > 2, see 
0) the bound identical to (|6.2|) also holds for the negative semi-axis. These bounds 
allow us to dehne the characteristic scales L^_i(F). Consider the case of L^_i{E) 
(the other one is analogous). By (|1.7|) , L^_i{E) is determined by the equality 


From 


we get 


||d>(n,F)|||,_^(^)=4||<h(l,F)-i|pT^ 


L+ i(F)3iog(--D < C'T^ 


and since ^(A) = A -|- 0(1) we can take 

~ , 61og(a;“^ 

L+_i(F) = OT^^^iW 
This proves the hrst part of Theorem 10. 


□ 










24 


ROWAN KILLIP, ALEXANDER KISELEV, AND YORAM LAST 


It remains to prove the second part of Theorem |1.6| , involving the lower bound on 
dynamics. The lower bound on dynamics for the Fibonacci Hamiltonian has been 
proved recently in [^. We present a sketch of the argument here, making explicit 
the behavior of the bound in the large coupling regime. The lower dynamical bound 
will follow from the continuity estimate for the spectrum. The idea is to study the 
behavior of solutions; as soon as appropriate bounds are available, one could apply 
the reasoning of Jitomirskaya-Last |^|, a recent result of Damanik, Killip and Lenz 
relating bounds on solutions to the a-continuity of the whole line operator, or 


Theorem 

Lemma 6.1. For every energy E in the spectrum of Hx, every solution u{n,E) of 
the equation {Hx — E)u{n, E) = 0 satisfies 

(6.3) ||u(n,E)|U>C'L", 

where k = log 20 independent of A. 

Proof. The proof follows closely Proposition 10 of [^. As shown by Siito |^, the 
Fibonacci potential obeys V{qn + 1) = V{1) for n > 3 and I < I < Qn- For n large 
and 1 < / < gn- 2 , we also get V (2g„ + /) = !/+ g „_2 + /) = !/(g „_2 + /) =V{1), 
so that 


(6.4) 


V(l) = V{qn + l) = V(2q„ + 1). 


By Lemma 1 of |^, we have for any 2x2 matrix B with detH = 1, 

max{|TrH|||HT||, > ^||T|| 


■il^lP 


for any 2-vector T. Therefore, 

(6.5) 

Let <F(m,/c) denote the transfer matrix that takes {u{m + l),u{m))^ to {u{k + 
l),u{k))'^. By (1^ ), for 1 < / < g „_2 we have <F(/,gn + 0 = + 1), 

and, moreover, Tr <F(/, + /) = Xn- Thus, (|6.5|) implies 


\n{qn + / + 1)P + \u{qn + /)p + \u{2qn + / + 1)^ + |u(2g„ + OP ^ 

^ -(!«(!+ 1)P +wot) 


4max(l, X. 


for 1 < / < qn- 2 , and 

\n{qn+i + / + 1)P + \u{qn+i + OP + \n{2qn-\-i + / + 1)^ + |M(2g„_|_i + /)P > 

for 1 < / < qn-i- Notice that for any E in the spectrum of Hx, for every n either \xr, 
or \xn+i\ is less than 2. By combining the above estimates, we easily deduce that 


\u{n,E)\\ 


Qn + 5 — 


\u{n,E)\\ 


qn 


( 6 . 6 ) 


4 
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for large n. The estimate (|6.3| ) follows directly from (|6.6|) . 

On the other hand, lochum and Testard m have shown that 


□ 


( log \/5 


log n 


The number d is dehned as a product d = ab'^, where 6 = 2c + 1, c = sup„|a;„|, 
and a = max(c, 2). For A large, one can check using (|5.4|) that c < A + 2. Hence, 
d < (A + 2)(2A + 5)^. Therefore we conclude that 


(6.7) 

where 


L 

E 

n=l 




C2(A) = f!^^(logA + 0(l)). 

log(o; 1) 


Recall that a measure fi is called a-continuous if it does not give weight to sets of 
zero a-dimensional Hausdorff measure. 

Proposition 6.2. The spectral measure pi of the operator Hx is a{\)-continuous for 
every A > 0, with a (A) satisfying 

2k 

(6.8) a(A) > ^ =0(logA)-^(l + 0((logA)-^)). 

« + C2(A) 


Proof. Given the estimates (|6.3|) and (|6.7|) , for the Fibonacci operator Hx the Propo¬ 
sition follows along the lines of Alternatively, one may use Theorem 1 of 


which shows for any discrete Schrodinger operator H on whole axis that if every 
solution of {H — E)u = 0 obeys 

< \\u{n,E)\\L<C2{E)H^ 

for E E then the spectrum of H is a-continuous in S with a = 2qi/{qi + q 2 )- □ 

Now we complete the proof of Theorem p..6| . 

Proof. The fact that a-continuity of the spectral measure implies lower bounds on 
dynamics is well-known. The original result is due to Guarneri 0,0 , with later con¬ 
tributions by many authors. Assume that the spectral measure pi^ of the Schrodinger 
operator H corresponding to the vector cf is a-continuous. Then by Theorem 4.2 
of the measure p,^ can be represented as a sum of mutually singular measures 
-|- /i‘^2, where is uniformly a-H51der continuous (/i'^i(/) < G|/|“ for any 
interval I) while can be made arbitrarily small. Moreover, it has been shown 

in proof of Theorem 6.1, that 


( 6 . 9 ) 


e-*'“<l>,\\ldt < CLa 
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(in dimension one). Multiplying (|6.9|) by a and integrating with respect to the kernel 
^^- 2 alT^ we get 


Jo 

In particular, given any > 0, we can choose 02 of sufficiently small norm and Cg so 
that for L(T) = CgT^, 


( 6 . 10 ) 


2 _ 

T 


g-2VT||g-imj,||2 


L{T) dt < g. 


Application of ( |6.10|) and Proposition completes the proof. 


□ 


7. Multidimensional problems 


Our purpose in this section is to show how our results for tridiagonal operators 
can be applied to more general problems. We are particularly interested in discrete 
Schrodinger operators of the form A + 1/ on £^(Z'^), dehned by 

(7.1) {{ A + V)i/j){n) = ^ i/j{m)+ V{n)-ilj{n), 

\n—m\=l 

where C M. 

Given any self adjoint operator iP on a separable Hilbert space Ti and a vec¬ 
tor 0 in the domain of 77, the cyclic subspace spanned by H and 0 is dehned by 
Tip = {/(T7)0 I / G C'oo(IP)}. Here Goo(M) is the set of continuous (complex valued) 
functions on R vanishing at inhnity and ^ denotes norm closure in 7i. Tip is an in¬ 
variant subspace for H and by the spectral theorem, the restricted operator H \ Tip 
is unitarily equivalent to multiplication by the coordinate parameter on the space 
L^(R, dyUp), where /ip is the spectral measure of 0 (and H). We assume that /ip is 
not supported on a hnite number of points and thus Tip is inhnite dimensional. The 
evolution given by (O is thus conhned to Tip (namely, 0(t) G Tip for any t) and so 
to study the time evolution, it suffices to consider the restricted operator H \ Tip. 

For large classes of problems, the moment vectors {Ti’^0}^Q are well dehned and 
the set of their hnite linear combinations is dense in Tip, namely. 


(7.2) 


^ span{i7"V'}^-o. 


We note that this holds for any 0 if T7 is bounded. If H is unbounded, then 0 must 
be chosen appropriately in order for this to occur (but such 0’s always exist, e.g., 
any 0 with a compactly supported spectral measure). We note that 0 must be in the 
domain of each of the moment operators and thus its spectral measure must have 
sufficient decay at inhnity to ensure hnite moments. This necessary condition is not 
sufficient, though, and determining whether this property holds or not for a given 0 
is in general a rather rich question, equivalent to the problem of whether the spectral 
measure of 0 can be determined from its moments. See for more information. 


The important thing for us here is that the property (Q) holds whenever H is of 
the form (|7.1|) on and 0 is any vector which is compactly supported on lA. 
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If the property (|7.2|) holds, then has a natural orthonormal basis which 

is obtained by applying the Gram-Schmidt orthonormalization procedure to the set 
of moment vectors Each has the form p„ = where Pn(-) is a 

polynomial of degree n. The polynomials {Pn(')}^o known as the orthonormal 
polynomials of the spectral measure The self-adjoint operator H \ is known 
P3| to have a tridiagonal matrix representation in the basis {pn}^o ^y choosing 
the phases of the p„ appropriately, this tridiagonal matrix can be made to have only 
real entries. Moreover, the off-diagonal entries cannot vanish. Thus, we see that 
can be viewed as being by identifying p„ with (5n+i, and then H \ has 

precisely the form (|1.2|) with a Dirichlet boundary condition. Theorems |1 . 1| , IL 


and are thus fully applicable in this general setting. We call the above tridiagonal 
representation of H \ the orthonormal polynomial representation. 

The problem, of course, is that when one is interested in wavepacket spreading 
dynamics, one is normally interested in how the spreading occurs in some natural 
coordinate space which is generally different from the orthonormal basis {pnjjJLo- 
The ability to apply our results to such a coordinate space depends on the ability to 
relate the vectors pn to it. Fortunately, for the case of an operator of the form ( [f.l|) on 
and a compactly supported initial vector ifj, an appropriate connection exists. 
For simplicity, take the initial vector to be the delta function vector at the origin 
of Z'^, namely V’ = <^o- Then one easily sees from (|7.1|) that H^Sq, and thus also pn, 
is supported in a ball of radius n -|- 1 in Z'^. Thus, any portion of the wavepacket 
which is supported on the vectors po, ...,pl must also be confined to a ball of radius 
L -|- 1 in Z'^. This means that any upper bound on the spreading of the wavepacket 
in the “coordinates” {pn}^o E immediately also bounding its spreading in Z'^. In 
particular. Theorem 0 takes the following form. 

Theorem 7.1. Let H be an operator of the form o on £^(Z'^). Let the character¬ 


istic scale L'f_i (E) be defined by ( [I.5|) for the orthonormal polynomial representation 
of H \ Hsq- Then for any T > 0 and L > 1, we have 

(7.3) {\\xL+ie-^^%f)T > Cpso ({E| L+_,(E) < L]) , 

where C is some universal positive constant and xl denotes the orthogonal projection 
on a ball of radius L around the origin in Z'^. 


We note that the scale [E) in Theorem [f.l| is dehned through the orthonormal 
polynomial representation, a fact which may make it seem like a somewhat obscured 
quantity in the multidimensional context. It is interesting to note, however, that 
even in the multidimensional context, it is actually among the more straight forward 
quantities to compute. Each of the moment vectors can be easily computed 

explicitly and it only depends on values of the potential within a ball of a correspond¬ 
ing radius n, and given the moment vectors up to some scale L, the corresponding 
Po,...,pi are obtained from them by the Gram-Schmidt procedure. The matrix el¬ 
ements of the (tridiagonal) orthonormal polynomial representation of H \ TLs^ are 
then easily obtained up to a corresponding scale L and from them one can explicitly 
compute solutions of equation ( |1.3|) up to a corresponding scale. It is thus possible, 
for any E and T, to compute (or more precisely, to bound within two integer values) 
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the scale by a straight forward computation consisting of hnitely many steps 

and only depending on values of the potential within a ball of a corresponding radius. 
This should be contrasted with the fact that multidimensional representations of gen¬ 
eralized eigenfunctions, such as the ones considered in the Kiselev-Last lower bound 
25| , depend on the entire potential and cannot be computed by a finite procedure. 


Acknowledgment. We would like to thank Svetlana Jitomirskaya and Barry Simon 
for useful discussions. 
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